Continuously Tunable Charge in Andreev Quantum Dots 
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We show that a quantum dot connected via tunnel barriers to superconducting leads traps a 
continuously tunable and hence fractional charge. The fractional charge on the island is due to 
particle-hole symmetry breaking and can be tuned via a gate potential acting on the dot or via 
changes in the phase difference across the island. We determine the groundstate, equilibrium, and 
excitation charges and show how to identify these quantities in an experiment. 
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The superconducting phase ip is usually associated 
with a dissipation- free flow of current; a phase twist by 
2tt within a bulk superconductor drives a magnetic vor- 
tex [l[ and a phase drop ip across a junction drives a 
Josephson current [2| . Only much later it has been rec- 
ognized that vortices carry a (small) charge as well 
an effect which is due to particle-hole symmetry break- 
ing. In this letter, we show that particle-hole asymme- 
try also generates a localized charge in a metallic quan- 
tum dot coupled to a superconducting loop, a so-called 
Andreev quantum dot see Fig. Ufa). This charge 
is continuously tunable (and hence fractional) through 
the gate potential V g or via the phase bias tp across the 
dot. It manifests itself for the (even parity) ground- and 
doubly-occupied excited states, while the (odd parity) 
singly-occupied excited state exhibits an integer charge, 
see Fig. He); as a corollary, the excitation charges are 
non-integer as well. The fractional charge can be ob- 
served via measurement of the associated telegraph noise 
signal due to the stochastic occupation of the Andreev 
states or through its dependence on the flux threading 
the loop. The latter provides a realization of a flux-to- 
charge converter allowing to make use of the device as a 
magnetic-flux detector. 

A tunable charge in a mesoscopic setup, a ring with 
a (nearby or interrupting) quantum dot, has been dis- 
cussed in Refs. 0, Qi the origin of this fractional charge 
is easily understood in terms of the extended nature of 
the wave function, with only a fraction localized on the 
dot. This contrasts with the fractional charge discussed 
here, where the 'missing' part is completely delocalized 
in the nearby superconducting condensate. Our charge 
then reminds about the fractional charge associated with 
excitations in a superconductor which has been discussed 
in the context of charge relaxation [f| (or branch im- 
balance @) in non-equilibrium superconductivity; there, 
too, the missing charge is 'dissolved' in the supercon- 
ducting condensate. In the following, we first discuss in 
detail the origin of the tunable charge which we then 
determine quantitatively. We discuss quantum fluctua- 
tions in the trapped charge as well as corrections due to 
weak Coulomb interactions, see Ref. S for a discussion 



of strong Coulomb effects and the Kondo effect in an An- 
dreev dot. We end with a discussion of telegraph noise 
allowing to detect the fractional charge. 

In superconducting mesoscopic devices the phe- 
nomenon of Andreev scattering [10( plays a central role: 
electrons (holes) incident from a normal metal (N) with 
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FIG. 1: (a) Andreev dot (SINIS) embedded in a supercon- 
ducting ring. Tuning of flux $ or gate potential V g serves to 
modify the charge on the dot. (b) The mixing of Andreev- 
and normal scattering at an SIN boundary leads to particle- 
hole symmetry breaking and to the appearance of a fractional 
charge on the dot. The phenomenon shows up in the vicin- 
ity of a resonance in the associated NININ device, here an 
electron type resonance, where the electron quasi-particle is 
multiply reflected while the hole component is not. (c) Charge 
trapped on normal- and Andreev dots versus position of reso- 
nance Se„; shown are the ground-state (gs) and single-particle 
excited-state charges for the singly and doubly occupied situ- 
ation, gs: as the resonance passes the Fermi level, a charge 2e 
is pushed out of the dot over a width F; the shaded areas in- 
dicate the importance of quantum fluctuations. Dashed lines 
(normal dot) indicate electron- and hole- type resonances; the 
associated trapped charges decay rapidly to the ground state 
charge (arrows). These resonances are transformed into eigen- 
states in the Andreev dot; they carry integer charge if occu- 
pied once (odd-parity) and a fractional charge when occupied 
twice (even-parity). 
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energies \E— e P | < A and scattering off a superconductor 
(S) are reflected as holes (electrons) with a Cooper-pair 
entering (leaving) the superconductor (here, e F and A 
denote the Fermi energy and the superconducting gap). 
In an SNS junction, electrons convert to holes at one 
NS boundary and back to electrons at the opposite NS 
contact, thus generating a particle-hole symmetric un- 
charged bound state (we ignore a small electron-hole 
symmetry breaking of order A/e F ). When the NS in- 
terfaces are augmented by thin insulating layers I in an 
SINIS junction, the particles in the normal-metal re- 
gion undergo imperfect Andreev scattering at the NIS 
boundaries with mixed Andreev- and normal scattering 
processes; the resulting Andreev levels involve unequal 
electron- and hole-like components and hence particle- 
hole symmetry is broken, see Fig. [TJb). Alternatively, 
the Andreev state can be understood as deriving from 
normal resonances in the associated NININ structure — 
replacing the normal leads by superconducting ones, a 
resonance residing in the gap region transforms into a 
sharp Andreev level with predominantly electron- (reso- 
nance above e F ) or hole character (resonance below e F ). 
The properties of these trapped Andreev states are easily 
tuned, either through shifts in their position within the 
gap via changing gate bias V g or via changes in the phase 
ip across the superconducting banks. 

The Andreev states give rise to new opportunities for 
tunable Josephson devices, e.g., the Josephson transistor 
111 . 1 1 21 ] ; here, we are interested in their charging proper- 
ties. Below, we consider small islands with a large sep- 
aration 6 n of resonances in the associated NININ prob- 
lem, 8 n > A, such that a single Andreev level is trapped 
within the gap region. Such Andreev dots have recently 
been fabricated by coupling carbon nanotubes to super- 
conducting banks [3, [3| , with the main focus on the 
competition between the Kondo effect and superconduc- 
tivity [13]. We are interested in sufficiently well isolated 
islands with a small width r„ of the associated NININ 
resonance, T n < A. Also, we wish to neglect charging 
effects in our first analysis here and hence demand that 
the Coulomb energy is small, Eq < T n . In summary, our 
device operates with energy scales Eq < T„ < A < S n , a 
situation which can be generated in the lab, see our es- 
timates below. In what follows, we determine the charge 
configuration of the ground- and excited states versus V g 
and if, find the thermal equilibrium charge and discuss 
the excitation dynamics of the Andreev states (due to 
elastic and Coulomb interactions) providing the charac- 
teristics of the telegraph noise in the charge state of the 
dot. 

The resonances in the NININ setup derive from the 
eigenvalue problem Ho^S — E"i> with Hq = —h 2 d 2 /2m + 
U(x)—e F and the potential U(x) = U pa (x+L/2) + U ps (x— 
L/2)] + eV g <d(L/2 — \x\)] describing two point-scatterers 
(with transmission and reflection coefficients \/Te ix and 
\f~Re %x ) and the effect of the gate potential V g , which we 



assume to be small as compared to the particle's energy E 
(measured from the band bottom in the leads), eV g <C E. 
Resonances then appear at energies E n — e^rm ~ x r ) 2 ; 
they are separated by 5 n = (E n+ i — E n _ 1 )/2 w 2E n /n 
and are characterized by the width r„ = TSn/nVR, 
where el = h 2 /2mL 2 . The bias V g shifts the resonances 
by eV g ; we denote the position of the n-th resonance rel- 
ative to e F by 5s n = E n + eV g — e F , see Fig. QIb). 

We go from a normal- to an Andreev dot by replacing 
the normal leads with superconducting ones. In order to 
include Andreev scattering in the SINIS setup we have 
to solve the Bogoliubov-de Gennes equations (we choose 
states with e > 0, see Fig. [TJb)) 



H A{x) 
A* (a) -H 



u n (x) 
v n (x) 



u n (x) 
v n (x) 



(1) 



with the pairing potential A(x) = A[6(-2 - L/2)e~ itp/2 
+@(x-L/2)e iv / 2 ); u n (x) and v n (x) are the electron- and 
hole-like components of the wave function. The discrete 
states trapped below the gap derive from the quantiza- 
tion condition (in Andreev approximation; the bar ' 7 ' and 
indices '±' refer to scattering probabilities and phases for 
the NININ setup and energies e F ± e) 



cos(S+ -S- -2a) =^JR+R- cos f3 + ^/T+T_ cosy. (2) 

The phase a = arccos(£*/A) is due to Andreev scattering 
and the phases S± = x±+k±L, k± — a/ 2m(e F ±e)/h ac- 
count for the propagation across the island; for symmet- 
ric barriers the phase (3 — (x+ — X+)~ (X- — X-) i s a mul- 
tiple of 7r and produces a smooth function y/R + R_ cos/? 
changing sign at every resonance 12, TB |. 

The solution of @ is straightforward and provides the 
spectrum shown in Fig. HJa). A small dot L < £ (£ 
the superconducting coherence length) results in a single 
Andreev state which lives near (but below) the gap A 
when the dot is tuned far off resonance with \Se n \ 3> A, 
e* s» A(l — r^/8fc 2 ); when the resonance resides in 
the gap region, \Se n \ < A, the associated Andreev level 
approaches the Fermi energy e F linearly, 



(i-r n /2A)|fe n |, 



(3) 



where we have assumed r„ <C A. As e* drops below the 
resonance width r„ the Andreev state becomes sensitive 
to the phase difference ip across the junction, 



4 « (1 - r„/2A) Vfe2+(lV2) 2 cos%/2). (4) 

The resonances above A follow closely the normal reso- 
nances e n , see Fig. [2fa). 

The Andreev state at carries a nontrivial 

charge which we find as the expectation value of 
the charge operator (for the space interval [—a, a]) 
Q a = e£ a dx'E (r=±1 iS>Ux)* t r(x) with = 

J2n[ U n( x ) a n,a + &Vn( x ) a n,-o] an ^ a n,<y, Oi\i a tne USUal 
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Bogoliubov operators; the J^n includes localized and 
extended states. Defining the quantities q® u = 
/" dxe\u n (x)\ 2 and q® v = J" dx e\v n {x)\ 2 (with the 
normalization q^ u + q^ v = e), we obtain the thermal 
equilibrium charge (we subtract the charge 2e(n—l) due 
to filled resonances and continuum states) 



<eq = 2/(4)<„ + 2[1 - /«)] (?»„; 



(•5) 



setting the Fermi function /(e) to zero we obtain the 
ground state charge g^ jgs = (0|Q a |0) = 2q® v , The exci- 
tation charge (of one spin state \a}) assumes the value 



<ex = HQ» - (0|Q°|0) = < u - ql, 



(6) 



These charges are localized around the dot, either strictly 
within the region [-L/2.L/2] as for the excitation charge 



a 



a ' 



or with an additional extension of size 



around the dot carrying a small weight [17| , q!^ eq « q^/eq 
(hence, the fractional charge is no£ due to a missing part 
of a wave function). The excitation charge can be di- 
rectly obtained from the voltage dispersion of the An- 
dreev state, q^x = ®V g s n'j combining this with the nor- 
malization condition, the other charges g£° eq and q^ gs 
follow immediately. The expressions ((3]) and (j4|) for the 
position provides the excitation charge, 



L/2 



e (1 - r„/2A)&„ 



Se n < r„, 



V<S4 + (r„/2) 2 cos2( v /2)' 

sign(fc„) e (1 - r„/2A), &„ < A, 
sign(<Je n )eAr£/4<$4) A < fe n . 



(7) 




FIG. 2: (a) Bound state (thick solid lines) and resonances 
e*' rcs (dashed lines) in the Andreev dot versus position Se n of 
the normal state resonance (tuned via the gate potential V g ); 
we choose F n = 0.2A. Also shown are the normal resonances 
e„ (thin solid line) and the dependence on the phase difference 
ip across the dot (inset), (b) Charges g^4x an d 9n,gB- Top 
inset: structure of telegraph noise with long periods of ground 
state (0) and singly occupied (1) states and the ground state 
trace interrupted by fast doubly excited states (2). 



The full numerical result is shown in Fig. [2jb): the 
charge rises linearly (with slope w 2e/r„ cos(y>/2); note 
the sharp change in charge for a 7r-tuned Andreev dot 
at ip — 7r) as the normal resonance crosses the Fermi 
level and then saturates at a value e (1 — r„/2A). As 
the normal resonance leaves the gap region, \Se n \ > A, 
the charge q^x vanishes oc AT^/Se^. We see, that tun- 
able fractional and localized ground state- and excitation 
charges appear on the Andreev dot each time a normal 
resonance crosses the Fermi level; this is the main re- 
sult of our paper. Note that the addition of the (frac- 
tional) ground-state- and excitation charges conspires to 
produce the integer charge of the singly-occupied odd- 
parity excited state, (a\Q L/2 \a) = q££ + q^'l w e. 

The above non-trivial charges are a consequence of the 
broken particle-hole symmetry and result from a super- 
position of states with definite integer charge, electron- 
and hole-like. As a result, the charges will undergo 
quantum fluctuations which are quantified by the vari- 
ance Sq a = [({Q a ) 2 ) - (Q a ) 2 } 1/2 . Here, we are in- 
terested in the charge fluctuations of the ground- and 
singly/doubly excited states |0), | ± 1), and |2); with 
(u\{Q a ) 2 \v) = £ M |HQ Q |^)| 2 the determination of Sq a 
involves the calculation of matrix elements (v\Q a \^i). A 
detector will measure the charge over some time inter- 
val r such that the relevant charge is given by the time 
average Q a = (dt / 'r)Q a (i); as a result, only ma- 
trix elements connecting states with energy difference 
less than H/t have to be considered; assuming typical 
measurement frequencies 1/r <C A/fi, we can restrict 
the sum over intermediate states to the four states 
|0), | ± 1), |2). The only relevant matrix element then 
is (0|Q a |2) = 2e f dxu n (x)v n (x) which is of order e. 
Hence, the charges of the singly-occupied states | ± 1) 
do not fluctuate, while the charges of the ground- and 
doubly-excited states fluctuate strongly, see Fig. [He). 

Next, let us comment on the effects of Coulomb in- 
teractions. Two important issues are the screening of 
the dot's charge and the mixing between charge states. 
Screening due to additional charge-flow to and from the 
Andreev dot is governed by the density of states. The 
relevant energy scale is given by the distance S n be- 
tween resonances (as these carry large local charges) and 
hence on-dot charge screening is irrelevant if Eq <C S n . 
Mixing of states is again governed by the matrix ele- 
ments (v\(Q a ) 2 \v)\ restricting the analysis to the four 
states |0), | ± 1), |2), we find that mixing occurs be- 
tween the ground- and doubly-excited states but is small 
if Ec <C e£. The energy scale Ec ~ e 2 /2C can be es- 
timated via the capacitance C w eL of the Andreev dot 
with e the dielectric constant; on the other hand, the dis- 
tance between resonances S n — hv F /2L depends linearly 
on L, too. The crucial dimensionless device parameter 
is the ratio 6 n /Ec = hv F e/e 2 ; assuming typical values 
e ~ 10 and v F ~ 10 6 m/s, we obtain a ratio 6 n /Ec ~ 30, 
allowing for a sequence of energies Eq < T n < A < S n as 
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required in our setup. In Ref. 1J| an Andreev dot in the 
form of a nanotube with L w 500 nm has been studied, 
that corresponds to a charging energy E c of order one 
Kelvin, hence the above inequalities can be satisfied in a 
realistic device. 

The fractional charge of excitations can be detected by 
various means; here, we discuss its observation through 
the measurement of the telegraph noise which arises due 
to the thermal occupation of states (we denote temper- 
ature by 9 and set k B — 1). The noise pattern in- 
volves the three states | ± 1) and |2); these states are 
filled and emptied through phonon absorption and emis- 
sion processes [l8]; in addition, fluctuations of the gate 
potential provide a competing channel. Two types of 
processes contribute to the telegraph noise, those cou- 
pling the ground- and singly-occupied states with rates 

701 and 710 = 701 exp(e/0), and those coupling the 
ground- and doubly-occupied states with rates 702 and 
720 = 702 exp(2e/0) (we introduce the shorthand e — £*; 
the rates connecting | ± 1) and |2) agree with the single- 
particle rates, 701 = 712 and 710 = 721)- The transitions 
1 involve a second quasiparticle in the continuum 
with energy E > A. 

The rates 701 and 702 arc derived from Fermi's Golden 
Rule with the perturbation given by the electron-phonon 
interaction H = g J dx n e (d x u) with u the displacement, 
n e the electronic density, and g the electron-phonon cou- 
pling assuming values of order 1 eV (we consider one- 
dimensional modes for both electrons and phonons); an 
alternative channel is provided by the fluctuating gate 
potential with H — J dxen e V g . Concentrating on elas- 
tic modes, the determination of 702 is straightforward 
and a simple estimate provides the result (N 2e is the 
Bose factor for the phonon state, a and k > 1/L de- 
note the lattice constant and the phonon wave num- 
ber) 702 ~ (g 2 /hmvl){a/kL 2 )N 2£ - As k < 1/L, we 
enter the zero-dimensional case and the result crosses 
over to 702 ~ (g 2 /frmv 2 ) ak N2 e . On the other hand, 
fluctuations of the gate potential contribute the rate 

702 ~ (e 2 /hC)N2 € - The numerical estimate with g ~ 1 
eV, mv 2 ~ 1 eV, and L w 500 nm provides the result 
702 ~ 10 12 s^ 1 N2 £ at k ~ which has to be com- 
pared with the Coulomb term 702 ~ 10 11 N2 £ s _1 . 

The calculation of the rate 701 involves an additional 
sum over continuum states which is dominated by en- 
ergies E ~ A; here, we estimate the rate for the si- 
multaneous occupation of the Andreev state and a con- 
tinuum state under phonon absorption and find 701 ~ 
{g 2 T/hmv 2 )(as/Lv F )y r F/Ae- A / 9 [l + e- £ /% where s de- 
notes the speed of sound. Note that kL ~ AL/hs ^> 1 
in the present case. Fluctuations of the gate potential 
contribute with a rate 701 ~ tT 1 {e 2 /C) ^6/Aer A l e [l + 
g-e/0]. usm g typical parameters 9 0.1 A sw 1 K and 
v F /s ~ 10 3 , we obtain numerical results of order 701 ~ 
10 10 s- 1 Texp(-A/6') and 7o i ~ lO 1 ^" 1 exp(-A/0) due 
to elastic and Coulomb interactions, respectively [lj|. 



Summarizing, we find that the processes <-> 2 are al- 
ways fast, while the <-> 1 transitions involve the expo- 
nential factor exp(— A/0), allowing us to reduce the rate 
701 dramatically by lowering the temperature. Hence, 
the <-> 1 processes can be effectively separated from 
the <-> 2 transitions, an effect which can be exploited 
in the experimental detection of the fractional charge. 

Given today's performance of single electron transis- 
tors (SET) operating with a sensitivity of ~ 10 _5 e/vTIz 
at frequencies / < 10 9 Hz [20j], the telegraph noise due 
to the thermal population of the Andreev state may be 
observable and thus the fractional charge on the Andreev 
dot may be measured. However, the processes <-> 2 may 
well be too fast, beyond today's time resolution of a SET. 
In this case we suggest two alternative schemes for the 
measurement of the excitation charge: i) Measuring the 
time averaged equilibrium charge g^ oq once at high tem- 
peratures, q£ jeq (0 » e) = q C n_ u + q% )V , and another time 
at low temperatures, q^^qifi <C e) = 2q® „, their differ- 
ence provides the excitation charge q% GX = q® u — q® v . ii) 
A measurement at low temperatures resolving the slow 
<-*• 1 transitions and averaging over the fast 0^2 
processes, see Fig. [21 measures the following (average) 
charges: given the occupation probabilities po = (1— f £ ) 2 , 
Pi = 2/ s (l — / e ), and P2 — f 2 , the average over the 
rapidly fluctuating regimes involving the states and 2 
provides the charge (g£)o,2 ~ 2q£ <v + 2/|^ ex , while the 
singly charged regime appears with the trivial average 
(?n)i = 1n,v + 1n,u- ^ l° w temperatures we have / £ «1 
and the difference - (g£) 0)2 « q%, u -q% >v gives q^ cx . 

In conclusion, we have discussed the peculiar prop- 
erties of the trapped charge on an Andreev dot. We 
have found continuously tunable fractional and fluctu- 
ating charges for the ground- and doubly-excited states, 
while the odd-parity singly-occupied state exhibits an in- 
teger charge. The transition dynamics of these states 
allows to observe these charges in the telegraph noise. 
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